ON THE STEADY STATES OF WEAKLY REVERSIBLE CHEMICAL 

REACTION NETWORKS 
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Abstract. A natural condition on the structure of the underlying chemical reaction network, 
namely weak reversibility, is shown to guarantee the existence of an equilibrium (steady state) in 
each positive stoichiometric compatibility class for the associated mass-action system. Furthermore, 
an index formula is given for the set of equilibria in a given stoichiometric compatibility class. 
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1. Introduction. The goal of chemical reaction network theory is to for- 
mulate conditions under which the dynamical fate of composition trajectories can 
be ascertained, even in complex reactions. Chemical reaction systems considered in 
practical applications are often very complicated. In principle, the number of species 
involved can be arbitrarily high. In practice, a full system can involve tens or hun- 
dreds of species. Although models with a handful of species are usually used, these 
systems are still "high" dimensional from the perspective of dynamical system. More- 
over, the reactions between the complexes is often only known approximately, since 
the determination of rate constants is quite difficult, if not impossible. A natural 
approach is therefore to attempt to give a qualitative description of the behavior of 
the reaction system. 

The modeling of chemical reactions can be achieved via the mass-action assump- 
tion. A typical chemical reaction network (CRN) obeying mass-action kinetics consists 
of three components: {S,y,K = (%)}, where 

1. S — {Si, S2, Sn} is the set of species involved in the chemical reaction. 

2. y = {Yi,Y2, ...jYm} is the set of complexes in the given CRN, with each 

g = l,2,...,m. 

3. If we have a reaction Yi Yj, then kij = ky^^Vj > is the rate constant for 
the reaction Yi ^ Yj, assuming mass-action kinetics. If there is no reaction 
from Yi to Y, , then we define kij = ky^^y. = 0. 

In this framework, the equation governing the evolution of chemical concentra- 
tions is given by: 



(1.1) x^f{x)=Y,h,x^^{Y,-Y,) 

where x g i?" = {x = (xi, . . . , Xn) '■ Xi > 0, i = 1, 2, . . . , n}, and each Xi is the con- 
centration of species Si in the chemical reaction network. 

The system we obtain is an array of ordinary differential equations with poly- 
nomial vector fields determining the behavior of the concentration of each species in 
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the reactor. This ODE could be quite complicated: the dimension of the equation 
could be dauntingly high, while the number of parameters (rate constants, for ex- 
ample) could be equally large. It is therefore quite surprising, but nevertheless true, 
that many system exhibit relatively simple dynamics. In fact, based on considerable 
experience, chemical experimentalists have developed an "intuition" that a "normal" 
chemical experiment will quickly lead to some steady state. 

So one of the primary concerns of chemical reaction network theory, up to now, 
has been to determine the capacity of a given reaction network system for multiple 
steady states. This theme emerged from the early work of chemical reaction network 
theory ([6], [8]), where chemical reaction networks were classified according to their 
deficiency. Horn, Jackson and Feinberg give a quite complete picture of the behavior 
of the deficiency-zero and deficiency-one mass-action systems ([1], [5]). Introductory 
material for chemical reaction network theory is given in ( [5] , [3] ) , and we will follow 
the notation and definitions in these papers as much as possible. 

The capacity of a reaction network for multiple states depends heavily on the 
underlying algebraic structure of the reaction system, given by the "react to" relation 
between the complexes. Two complexes are linked together if there is a reaction 
between them, and we say two complexes are in the same linkage class if we can find 
a finite sequence of reactions between them. It is not hard to check that the linear 
subspace spanned by all the reaction vectors Yj — Yi, for each reaction Yi — > Yj, is 
invariant under the flow induced by (|l.ip . It is called the stoichiometric subspace. 
Each linear manifold parallel to it is called a stoichiometric compatibility class{SCC)^ 
which is also invariant under the flow. 

The "react to" relation between the complexes also gives us various notion of 
reversibility, as alluded to above. If each reaction A ^ B in the chemical reaction 
network is accompanied by the reverse reaction B A, then this reaction network is 
called reversible. If for each reaction A — i? we have a reaction chain from B to A, i.e., 
we can find a finite sequence of reactions B — Ci, Ci — C2, . . . , C„_i — )■ C„, C„ — )■ A, 
then this reaction system is called weakly reversible. We would consider solely weakly 
reversible chemical reaction networks in this paper. 

In the early eighties, Nachman obtained the existence of steady states in each 
stoichiometric compatibility class if there is only one linkage class in the chemical 
reaction network ([ID]). Using this result, Feinberg showed the existence of steady 
states in each SCO (personal communications) if the chemical reaction network is 
weakly reversible and the deficiency is equal to the summation of the deficiencies of 
each linkage class. So the next natural question is about the existence of steady states 
for general weakly reversible systems. 

We are interested in the following two problems: 

1. Does there exist a positive steady state for (|l.ip ? 

2. Does there exist a steady state in each nonempty, positive stoichiometric 
compatibility class for 

The main result in this paper gives a affirmative answer to the two questions 
above, under a natural condition. 

Theorem 1.1. For each weakly reversible chemical reaction network obeying 
mass-action kinetics, the flow of has finitely many (at least one) steady states 

in each nonempty, positive stoichiometric compatibility class. 

Our result shows that the experimental chemists' "intuition" is well grounded, 
at least for weakly reversible and mass-action chemical reaction networks. It shows 
that for such systems there will always be some steady state in each stoichiometric 
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compatibility class(SCC), and although the numbers of steady state may differ for 
different SCC, the summation of the indices for the steady states remains the same 
(see Corollary 5.2). Moreover, this does not depend on the specific values of the 
positive rate constants, as long as the system remains weakly reversible. This would 
be quite appealing to experimental chemists since the rate constants are difficult to 
measure. 

Remark 1.2. The component of the complexes Yi is called the stoichiometric 
coefficients, and for typical chemical reactions, it should be a nonnegative integer. 
But our result holds when the stoichiometric coefficients are real, thus for example 
we allow the following reaction: 

0^2A + B, 2.5B + C^0, -2C + D = E 

as long as the chemical reaction is weakly reversible. Therefore the case where there 
is "source" and "sink" for the chemical reaction system can still be described by our 
result. 

Plan of the paper. In section 2, we define two matrices C, R that determine the 
reaction network structure, discuss their properties and reformulate the existence 
problem as that of an intersection problem of two hypersurfaces. In section 3, we 
restrict ourselves to the case of one linkage class. A vector-valued function G{z) is 
defined, and a priori estimate of (G(z), z) is given and used to solve the intersection 
problem in a special case. In section 4, we discuss the case of I linkage classes and 
construct a bounded, convex set so that the intersection problem is reduced to proving 
that a corresponding vector field is pointing in at each point of its boundary. Thus we 
give the proof of the existence of steady states in the whole i?" space. In section 5, we 
utilize the remark given at the end of section 4 to prove the existence of a steady state 
in each stoichiometric compatibility class, and then we use a homotopy argument to 
give the proof of the index formula. 

2. Reaction netvifork structure and reformulation of the equation. It 

has long been known that the existence of steady states depends heavily on the al- 
gebraic structure of the chemical reaction network. For example, when the reaction 
network is reversible, forest-like and the deficiency of the network is 0, one can show 
the existence of steady states in i?" (the reader is referred to ([9 ) concerning the 
precise exposition). But it is not obvious how to extend the result to the weakly 
reversible reaction networks of arbitrary deficiency. 

The information about a weakly reversible reaction network can essentially be 
decomposed into two parts: one about the configuration of the complexes in the i?" 
space, offered by the set y = {Yi, I2, • • ■ , ^m}; the other one about the reaction 
information between the complexes, given by the rate constants kij, i,j = 1, 2, . . . , m. 
We give the following 

Definition 2.1. For each weakly reversible chemical reaction network {5, 3^, K}, 
the configuration matrix C is given by 



C = 



( 


] 






2/12 • 


■ Vln \ 


Y2 






2/21 


2/22 • 


■ V2n 


\ Ym 


) 


7nxn 


\ J/ml 


J/m2 • 
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while the relation matrix R is given by 



R. — {Tij)rnxm! witk Tij — 



kij , if i ^ j. 

^ 'l2s^i ^is, if i — j- 
Therefore the reaction network structure is completely determined by the two 
matrices {C, R}. Notice that the property of weak reversibility is uniquely determined 
by the relation matrix R. They will both play a role in the existence of the positive 
equilibrium, but it will turn out that the specific structure of R will be crucial in the 
following development. 

System (|l.ip can be rewritten in terms of {C, R}: 



(2.1) 



^In X-Y2 



) 1 X m 



( \ 

Y2 



R C 



where {C, R} are given as above. 



We list some basic properties of the relation matrix R. 



1. R = (r, 



I'j Jmxm 



satisfies: 



Z^jLi = 0> foi' « = l,2,...,m 
nj > 0, for i ^ j 
ru < 0, for i = 1, 2, . . . , m 

2. Every two complexes in the same linkage class are linked by "reaction chains" . 
Mathematically, this means for each pair Yi,Yj in a certain linkage class, 
there exist Yi^,Yi^, ■ ■ ■ ,Yi^ such that Yi^ — Yi,Yi^ — Yj and ri^_^i^ > for 
k — 2 , 3, . . . , s. 

3. If the CRN has I linkage classes, / > 1, then by relabeling Yi, 1^2, • ■ • , Y„i we 
can write R as the following diagonal block form 



(2.2) 



R 



/ Ri 



V 



R2 



Ri / 



with each R^ being the relation matrix for a linkage class. 
The only thing in item 1 which needs explanation is ru < 0. If r^; 



for some i, 



then rij = for j ^ i. Thus the complex Yi does not react to any other complexes. By 
weak reversibility of the chemical reaction network, no other complex will react to Yi 
either. Thus Yi can safely be discarded and we only need to consider the subnetwork 
without Yi. Therefore we will always assume that each complex Yi appears at least 
once in the chemical reaction network. 

Remark 2.2. It is interesting to observe that the classical Perron- Frobenius 
theory for matrices with positive entries can be used to get information about the 
eigenvalue distribution of the relation matrix R. For example, we can prove that R 
has only nonpositive eigenvalues, and the algebraic multiplicity of the zero eigenvalue 
of R is exactly the number of linkage classes within the CRN. Also, the eigenvector 
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corresponding to the zero eigenvalue is nonnegative. This spectral information seems 
to have decisive influence on the stability property of the phase portrait of the flow. 
The interested reader is referred to ([7]). 

Notation 1. For each CRN we define the norm of its relation matrix R to he 
|R| = maxi=i^2,...,m \rii\- Also, we define r(R) = niini-^j i ,,g{i_2^. ..,,„} {rij : > 0}. 

Our next task in this section is to transfer the existence problem to one of inter- 
section of two hypersurfaces. The idea is to see that x S i?" is an equilibrium point 
for system (|2.ip if and only if 

/ glnx-Yi glnx-Ya ^lnx-Y^\ ■ Ti G { Z G R"" : Z ■ C ^ 0} , 

V / 1 X m ^ J ' 

So it is natural to consider the intersection 

(2.3) {( e'"^-^! e^"^-^^ ... gi"^-^" ) • R : a; e i?" } n {z G i?" : ^ • C = 0} 

if we want to find positive steady states for (|2.ip . 

Notation 2. N = {z e : z ■ C ^ 0}, K = {z e R"" : z = x ■ for some x e i?"}. 

Note that K = N^, = N where and mean the orthogonal subspaces 
of N and K in i?™ respectively. 

Definition 2.3. We define the relation function G{z) : R™ — i?™ as 

(2.4) G{z) =^ • R = (e^^e^^ . . . z^")R. 
Then we have the following 

Proposition 2.4. There exists a positive equilibrium point for \2.1]] in i?" if 
and only if 

(2.5) G{K)nK^ ^d). 
Proof. : It suffices to notice that 

G{K) = |( e'""'-^i e'"^ '^^ 

and 

K-^ = N = {z e R"" : z -C = 0} 

then from equation (j2.3p and the preceding discussion, we obtain the result .□ 

Remark 2.5. The idea of using the intersection of two hypersurfaces to prove an 
existence result can be traced to Felix E. Browder in (J^). Although it turns out that 
the relation function G{z) : i?™ i?™ is not monotone, the fact that we are dealing 
with a finite dimensional ODE still pulls us through. 



ginx-Y^ \ .R : a; € Rl] 

I Ixm + J 
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3. The case of one linkage class. In this section we will solve the in- 
tersection problem, assuming that the CRN has only one linkage class. The general 
strategy is as follows: to prove that G{K) H ^ 0, it is equivalent to showing that 

e Hk o G{K), (*) 

where Ili^- : i?™ K is, the standard projection operator. One sufficient condition for 
(*) to hold is that there exists a large ball B{r) = {z € K : \z\ < r} K such that 

(n;^oG(z),z) <0 

for z G dB{r) = {z e A' : |z| = r}. Then by the Brouwer Fixed Point theorem we 
obtain (*). 

So basically, we need to estimate 

{T\koG{z),z) = {G{z),z) 
for z & K and try to show that we can find a ball B{r) such that 

(G(z),z) < 

for z S dB[r). But here we encounter some technical difficulties: one quickly 
finds that (G(z),z) will be identically zero if z lies in the 1 — d subspace D = 
{z = a • (1. 1, . . . , 1) € -R™ : a € R}. Thus we will have to avoid this degenerate sub- 
space D. 

Definition 3.1. We setl^ = (1, 1, . . . , l)ixm,Pm = {z e i?™ ■ km ■ z = 0}, and 
we define the norm on i?™ as jzjj = J2iLi Zi' ■ If z € Pm, we define another equivalent 
norm \z\ = maxi=i^2,...,m ^j- 

We start from the following 

Lemma 3.2. Let s > 1 be fixed. Then 

(3.1) max yi -I- 6^^2/2 + ■ • ■ + e^=-iys + e^=^(-£) ^ -oo 

— L<yi<0,i— l,2,...,s 

as L ^ +O0. 

Proof. We will show p.ip by induction on s. 
Step 1: If s =1, then 

max yi + e'^^{—L) = max {— L, —L + e~^{-L), -1 — Ini} — > — oo 

-L<yi<0 

as i — > +00. 

Thus the lemma is true for s = 1. 

Step 2: Now suppose that the lemma is true for s < k. Define 

Ak,L^ max {yi + ey'y2 + ... + €^"{-1)} 
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then by assumption we know that A^^l —oo as L — )• +00. 
For s = A; + 1, let —L < j/i, 2/2, • • • , < be fixed, then 

*(yfe+i;yi,y2,. . ■,yk) =yi+ey'y2 + ■■■ + ey-yk+i + ey^+^-L) 
will achieve its maximum at yu+i = 0, — L or j/^ — InL. So we have 



(3.%+i,L = max *(yfe+i;yi,y2,---,2/fe) < max{-L, Afe,L,Afc,i„L}, 

— Z/<yi<0,i=l,2,...,fe+l 

therefore Ak^ix —00 as L — > +00. Thus the lemma is true for s = fc + 1. 
Step 3: Combining Steps 1 and 2 we see that the lemma is true for all s > 1. □ 

The fundamental lemma is the following 

Lemma 3.3. If the CRN has exactly one linkage class, then there exists a function 
c(r) on with c(r) +00 as r ^ +00 such that 



(3.3) {G{z),z)<-c{\z\)e 

for all Z <E Pra- 

Proof. For all z S P^, we have 



1^1 



= e^-l^limR^;* = e^''-l^limR(^;* - |^| 1^) = e^'Rj/* 



(3.4) I 



where y = z — \z\ . 

Notation 3. Bl = {y G R"^ : -L < < 0,V« = 1, 2, . . . ,TO;minj yi = -L,maxj y 

Then lemma 3.3 is equivalent to the following 
Claim: Define F : Bl ^ R as F{y) = e^Ry*, y G Bl, then 

Cl = niax F{y) —00 

as L ^ +00. 

Proof of Claim: Given y G Bl, by relabeling Fi, I2, • • • , Ym we can write y as 

( -L<y,< 0,V1 <i<t2 
^ \ti+t2 + t3 = m,ti>l,t3>l ■ 
So we have 



F{y) = eVRy' = ( e^'i e^*^ e'^ia ) 
(3.5) 

< e^Ri2y* + e«R23(-il*) + e^Ri3(-il*) + e-^iR33(-il*) + 6^22^* 
^ , ' ^ « ' ^ , ^ 

(1) (2) (3) 



Rii 


R12 




of 


R21 


R22 




1 


R31 


R32 
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Estimate of term (3): we have 

r term{3) = ^ {- L)lt^R33ll + eiRaaf < h |R| Le'^ + |R| (max_oo<.<o \e's\) 
I <(t2+<3)|R|e-i<f |R| 

(3.6) 

To estimate terms (1) and (2) we need to consider two cases: 

• Ri3 ^ 0, by which we mean there is at least one element € R13 which is 
larger than 0, then 

term{2) < {-L)aij < -L ■ t(R). 

Notice that term (1) will always be nonpositive since each element in R12, R23 
is nonnegative, and each element in i/^ and —Ll^^ is strictly negative. 
So we have 

777 

(3.7) ierm(l) + (2) + (3) < — |R| - L • t(R) 

e 

• If R13 = 0, notice that 

(Rii R12 R13 \ 
R21 R22 R23 corresponds to the decomposition of the CRN to 
R31 R32 R33 / 

three subnetworks Ni = {Yi, 1^2, . . . , J , A^2 = {Yt,+i,Yt,+2, . . . , Yt,+t,} , N3 = 

= {Yt^+t2+ijYt-^+t2+2, ■ ■ ■ ,Ym}- Rii corresponds to the relation matrix for 
Ni,i = 1,2,3, (Although Ni together with R^^ may not correspond to a 
weakly reversible chemical reaction network) while Ry , i =/= j G {1,2, 3} gives 
the information about reaction vectors from Ni to Nj. 

Now Ri3 = means that each Yi belonging to A'^i must pass N2 to get to Yj 
which belongs to N3. So we can always find a "reaction chain" 

Y,,^Y,,^...^ Y,^,s> 3 

such that Yi^ € Ni, Yi^ e A^3, and Yi^ € N2 ioi k e {2, . . . , s - 1}, which 
means that ai^i^ j 01213 ) • ■ • : > with 

aiii2 G R12, ai2i3 7 O^isiiJ • ■ • J ais_2is-l ^ R22, G R23- 

So now we have 

term(l) + (2) + (3) < a,,,,y„ + ey^2a,^,^y,^ +... + e^'-i a,_,,^ (-L) + f |R| 

< r(R)[y,, + e^'2y,3 + gV-ay,^ + . . . + e^-^y,,_, + e^'-^ (-L)] + f |R| 

< r(R) • ^_2,L + T ' 
(3.8) 

where in the last inequality of p.Sp we have used lemma 3.2. 
Combining p.7p . p.8p we see that 
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— |R| - L • t(R), t(R) • As-2.L + — |R| L 
„ - „ e e J 

thus Cl — maxjjgs^ P{y) ^ ~oo as i — > +00. The proof of the Claim is now 

complete, and the proof of Lemma 3.3 follows from that of the Claim. 

□ 

Remark 3.4. This lemma was obtained by Adrian Nachman, around 1983. 
Using this lemma, Nachman was able to show the existence of positive equilibrium 
points for weakly reversible chemical reaction networks with one linkage class. After 
this Feinberg showed the existence of equilibrium points when the CRN is weakly 
reversible and has the property that the deficiency of the CRN is equal to the sum of 
the deficiencies of each linkage class subnetwork. 

Corollary 3.5. If G{z) = • R where R is the relation matrix for some CRN 
of only one linkage class, then for all C > 0, there exists r{C , R) > such that for 
all w € Pm, \w\2 < C, we have 

{G{z + w),z) <0, z e P„, |z|2 > r(C',R). 
Proof. First observe that 

iG{z + w),z) = 6^+"" • R • z* = • R{w) ■ z* 

where R(w) Diag{e'"\e'^\ . . . , e'"") • R, with Diag{e'"\e'''\ . . . , e""'™) being the 
m X m diagonal matrix. Thus R(w) is still a relation matrix of some CRN. 

Think about R(w) — Diag{e™^ ,e™^ , . . . ,e'^'") • R as a perturbation of R, in 

the class of relation matrix. Fix C > 0, then for w e |zz; S Pm,\w\2 < cj, the 

proof of Lemma 3.3 still holds for R(w), specifically the estimate p.9p holds for 
w e Pm, \w\2 < C, with R replaced by R(it;). Now observe that 

min t(R(w)) > 0, 

\w\^<C 

then due to the compactness of |u) e P„j, {wl^ < we can choose r{C\ R) > large 
enough such that 

iG{z + w),z) = (e^ • RH, z) < 

holds for w e G Pm, \w\2 — ^} uniformly, where z G Pm, \z\2 > r(C',R). □ 

Corollary 3.6. If G{z) = • R where R is the relation matrix for some CRN 
of only one linkage class, then there exists constant to (R) > such that 

(G(z),z) <to(R), zePm 

Proof. From Lemma 3.3 we know that there exists rp > 0, such that for all 
z € Pm, \z\ > To, we have (G(z),z) < 0. Let to — max(0, niax2£p_^.|2|<ro |(G(z),z)|). 
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Since (G(z), z) is continuous for z £ Bp^{rQ) — {z £ P„ : \z\ < rg}, we have < to < 
+00. □ 

From Lemma 3.3 we obtain easily 

Lemma 3.7. // the CRN has only one linkage class, then for any subspace H of 
Pm, we have that 

(3.10) G{H)r\H^^% 

Proof. Since R™" = H , we can define the projection operator Hh : R™ H 

as 



(3.11) TIh{z) = zi 

where z = zi + Zi, z\ £ H,Z2 & . 

Then G{H) n iJ-^ 7^ if and only if e IiH{G{H)). From we know that for 
r sufRciently large, 



(Hh o G{z), z) = (G(z), z) < -c(|z|)el"l < 

for z e Snir) = {z e H : \z\ = r}. 

By Brouwer's Fixed Point Theorem, there exists zq € Bnir) = {z £ H : \z\ < r} 
such that Hh o G{zq) = 0. Thus we have 



e nH{G{H)) 

which means G{H) n i/^ 7^ 0. □ 

4. The case of I linkage classes. In this section we will discuss the case of / 
linkage classes. Without foss of generahty we assume the relation matrix R has block 
form (j2.2p . and each is a x m.i matrix for i ^ 1,2, ... ,1. 

The content of this section is divided into two parts: the first part is about 
the estimate of {G{z),z). Similar to the case of one linkage class, there exists an I 
dimensional degenerate subspace Di for G, i.e, {G{z),z) = restricted to Di. We will 
give an estimate of (G(z), z) restricted to the orthogonal subspace P„i,i of Di (Lemma 
4.3), and correspondingly show that G{H) n ^ when H C Pm,i (Lemma 4.4). 

DEFINITION 4.1. Let Pm,l ^{ze i?™ : = 0, E'ZVZ = • ■ • > J:T=rn-m, + l = 

Since Pm,i C Pm, \z\ is still a norm for z G Pm,i- 
For z — {z^^\z^'^\ . . . , z^')) e Pm.i where 

^ — (^1 1 ^2 7 ■ ■ • ; ^mi ) 1 ^ — (^mi + 1 7 ^mi+2 i ■ ■ • 7 ^mi +m2 ) 7 • ■ • 7 (^m— m; +1 7 ^m— mi+2 7 • • • 7 

we define Hi : Pm.i ^ Pnii as 



n,(z) = z«, z = l,2,...,L 
Remark 4.2. By Lemma 3.3, there exists < Rq < +oo such that 
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(4.1) 



for all G P„ 



> Ro,i = 1,2, 



We have the fohowing 

Lemma 4.3. If the CRN has exactly I linkage classes, then there exists a function 
c*{r) on with c*(r) — > +oo as r — > +cxi such that 



(4.2) 



(G(z),z)<-c*(|z|)el' 



for all z G Pm.i- 

Proof. Using the estimate in Lemma 3.3 for one hnkage class we have 



(G(z),z) = (e^'" e^"' ... e^'")R 



I + e^"R,z«* < _ min{ci(|z|), C2{\z\), . . . , q(|z|)} el^l + ^Li ^» < -c*{\z\)e\^^ 
(4.3) 

where c*(r) = min {ci(r), C2(r), . . . , Q(r)} — ; with Ci{r),ti corresponding to 

the relation matrix R^, ? = 1, 2, . . . , Z, respectively. □ 
From Lemma 4.3 we obtain easily 

Lemma 4.4. // the CRN has exactly I linkage classes, then for any subspace H 
of Pm,i, we have that 



(4.4) G{H) r\H^ 

Proof. Since i?™ = H®H^, we can define the projection operator Hh ■ -R'" H 



as 



(4.5) 



IIh{z) = zi 



where z = zi + Z2, z\ G Pl,Z2 G i?^. 

It follows that G(B) n 7^ iff G ^h{G{E)). From (gSJ we know that for 
r sufficiently large, 

(n^f o G(z), z) = (G(z), z) < -c*(|z|)el^l < 

for z G Snir) = {z £ H : |z| = r}. 

By the Brouwer Fixed Point Theorem, there exists zq G BH{r) = {z E H : \z\ < r} 
such that Hh ° G{zq) = 0. Thus we have 

OellHiGiH)) 



which means G{H)nH^ ^9. D 
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Now we turn to the second part of this section. We want to discuss an arbitrary 
subspace K of i?™ and show that (|2.5|) is true. The natural idea is to transfer the 
intersection problem for K to that for Pm,i- Therefore the possible configuration of K 
in i?"' = Di(BPm,i is first discussed, then via Lemma 4.5 the intersection problem (j2.5p 
is transformed to (|4.6I) . which is an intersection problem with respect to Pm.i- This 
is achieved at the expense of adding an arbitrary linear map F from some subspace 
of Pm,i to Di, which leads to the change of relation function G to that of G*. 

But the essence of this section is actually Lemma 4.5, so a few words about 
the idea of proof are helpful. The idea is to construct a domain ilnir) C H which 
resembles the ball Bnir) with a slight modification. We want to choose VLuir) so that 
it is compact and convex and thus homeomorphic to the ball. Moreover, we want to 
show that on the boundary of VLuir) we have 



where z G dflnir), n{z) means the outward normal of dnnir) at z. 

Then the Brouwer Fixed Point Theorem can be applied to obtain Lemma 4.5. 

Lemma 4.5. For any subspace H C Pm,i, o,nd for any linear map F : H ^ Di, 
we have 



or, equivalently, there exists x £ H, such that G{x + F{x)) E . 

Proof. For the sake of simplicity, we give the proof for 1 = 2. The proof for I > 3 
case is similar. 

When I = 2, we define 

(4.7) Hi=Hnn^\o), H2 = Hnn^^o), H3 = h n {Hi + H2)^ 

where Hi : Pm,i — > Prm , « = 1, 2 is as given in Definition 4.1. Then H = Hi® H2® H^ 
and H,±Hj for i^je {1, 2, 3}. 
Thus given z H, we have 



where zi = (0,xi) € Hi, zi = (a;2,0) G H2, Z3 = (2/1,1/2) G ^^3- Then there exists 
Ci > C2 > such that 



for all Z3 = {yi, j/2) G H3, since by construction Hi restricted to H3 are isomorphisms 
for i = 1,2. 



{Uh o G*{z),h{z)) = {G*{z),n{z)) < 



(4.6) 



Go{I + F){H)r\H^ ^% 



(4.8) 



Z = Zi + Z2 + Z3 = (0, Xi) + {X2,t)) + (2/1,2/2) 



C2I2/2I2 < \Vl\2<Cl\V2\2 
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First, we have the fohowing 

Claim 4.6. If r > C^^\ then we have that flHir) is homeomorphic to Bh{1) — 
{z (z H : \z\2 < 1}; thus dflnir) is homeomorphic toSni^) ^ dBni^) = {z ^ H :\z\ — 

Proof of Claim 4.6: ^h{t) is homeomorphic to -B/f (1) due to the fact that for r > C*^^^ , 
flnif) C is a compact and convex set, with nonempty interior, thus is homeomor- 
phic to Bh{1)- The second claim follows from the fact that the restriction of a 
homeomorphism is also a homeomorphism. Q.E.D. 

Notation 4. We define the structure function G* : H ^ P^.i as 

(4.9) G*{x)^G{x + F{x)). 

For the linear map F : H ^ Di, we will also write F as 

F{z)^{F, (z) 1„^ , F2 (z) 1„^ ,...,Fi (z) 1„, ) , 

where Fi(z) : H ^ R is a linear function for i = 1,2, ... ,1. 
Next we show that 

Claiim 4.7. For r > sufficiently large, we have (IIh o G*{z),h{z)) < 0, for all 
z € 9ri^f(r), where n{z) means the outer normal offlniT) at z d c?r2/f(r). 
Proof of Claim 4.7: When r > is sufficiently large, we have d^lnir) = S'iUS'2U53US'4, 
' S^ = LeH:\z\,=r,\z,\,<r~^,\z,\,<r~^,\z,\,<r~^\, 

S2^{zeH:\z,\,^r-^, \z\, < r, jz^l^ <r~^, \z,\, < r - ^} , 
where < > (2) (2) (2) i 

Ss^lzeH: IZ2I2 = r- ^,|zi|2 < r~ ^,|z|2 <r,|z3|2 <r-^|, 

^ S, = {zeH:\z,\, = r-^,\z,\,<r~^,\z,\,<r~^,\z\,<r} 
with outer normal n(z) = tA-, t^h-, t^Vi respectively. 
Thus for z G S'l, we have 

(Hh o G*{z),h{z)) = (Hh o G*{z), j^) - (G*(z), j^) 

Fl^ I^l2 

Now from z G S'l we have |z|2 — r, \z1\2 <r~ \z2\2 l^f— ^v-, thus 



U l/^ \^^\2 + k3l2^ = \z\2' - kil/ >r^-{r- ^? = 2C(^) C(^) 
^ ■ r l-il/ + l-3l2^ = \z\2' ~ >r-~{r^ = 2C(^) - ^ > C(^) 

for r > C(2). 

• If k3l2 < then from (|4.1ip we have 



\X2 + yi|2' > |2;2|2' = k2|' > ^ > i?o' 



(4.12) 

since H1J.H3, H2-LH3. Thus by Remark 4.2 and (|4.10p we have 
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• If [zalj^ > then we have 



(4.14) 



Thus, similarly, we also have that (|4.13p holds. 
Summarizing the estimate above, we have shown that if r > C'^^\ then 



(4.15) 



for aU z G Si. 

For 2 e 5*2, we have 



{nHoG*{z),n{z)) <0 



(4.16) 



{Hh o G*{z),n{z)) = [Wh o G*(z), ^) - (G*(^), ^) 

|zil2 N1I2 



Now, from z € 52 we have |zi|2 = r — ^7—, l^^lj < thus we have 



(4.17) 



\y2\,'<\z^\l<\z\l-\zi\l<r^-{r-^f 



which means \y2\2 < v2C^. 

With Corollary 3.5 of Section 3 in mind, we choose ri > satisfying ri — > 

r(V2C(2), R2)), such that when r > ri, we have \x1\2 — \zi\^ — r — > — 
thus combining the two inequalities above we have 



e-^+^^+^^(-)1^2R,(xi)* ^ e^^(-)-e-^+^^R2(xi)* ^ ^ 

kll2 I^ll2 



Thus we have shown that if r > ri , then 



(4.19) (nHoG*(z),n(z)) < 

for all z G 82- 

Similarly for z G 5*3, we can choose r2 > satisfying r2 — > r{V2C^, Ri)), 
such that when r > r2, we have 
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e 



X2+Vi+Fi{z)l 



for all z £ S3. 

Now for z e 5*4, we have 



(n^ o G*{z),hiz)) = (n^ o G*iz), ^) = {G*iz), ^) 

\^3\2 I^3l2 



z e 5*4 implies \z3\2 = r — [zL < r, thus 



(4.22j 



Ixil^ = \z,\l <\z\l- IZ3I2' < - (r - = 2C(2) - ^ < 2C(2), 

I = |.2|^ < |.|^ |.3l2^ <r^-{r- = 2C(^) - ^ < 2C(^) 



which means that |a;i|2 < V2C(2), IX2I2 < V2C(2). 



By Corollary 3.5 of Section 3, we can choose ra > satisfying — > 



max(Wi^, v^l + Ci^) •max(r(V2C(2),Ri),r(V2C(2),R2)), such that when r > 



we have IZ3I2 = r — > — which leads to 



(4.23) 



2;il2> JTS^k3l2>K\/2CW,Ri) 



^ 1^212 >y^T^k3|2>KV2CM,R2) 

Thus, we have that when r > r^, 



(4.24) 



{Hh o g*{z), n(z)) = :z^2(hl 

_ e^i<').e^2+«lR^fa^)t+e^2(»).e^i+«2R3(j^3)t 

^3l2 



for all z G 5*4. 

So combining (|4.15p . (|4.19p . (|4.20p and (I4.24p we have finished the proof of Claim 
4.7. Q.E.D. 

Combining Claim 4.6 and Claim 4.7, also noticing the discussion before Lemma 
4.5, we finish the proof of Lemma 4.5. 
□ 

Remark 4.8. It is easy to see that Hh o G* is an analytic function defined on 
H, so there can only be finitely many z S Oij(r) C H such that Hh o G*(z) = 0. 
Also from the proof of Lemma 4.5 we can see that outside flnir) there is no z £ H 
belonging to G*{H) n H^. 



Remark 4.9. It is an interesting observation that in the proof of Lemma 4.5, we 
never use the fact that F : H Di is linear. In fact, the reader may readily check 
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that if we change the word "linear" to "nonlinear" in Lemma 4.5, the proof of Lemma 
4.5 still goes through. This observation will play a fundamental role in the proof of 
Theorem LI. We write is as the following 

Corollary 4.10. For all H C Pm,i and for all continuous map F : H ^ Di, 
we have 



(4.25) G*{H)nH^ = Go{L + F){H)nH^ ^(D 

or, equivalently, G{x + F{x)) G for some x E LI . 

Lemma 4.11. For all H' C i?" such that LI' n Di {6}, we have 



G{H') n i7' ^ 0. 

Proof Let H = LEp^ , {LI'). From the assumption LI' n Di = H' (1 , = {9} we 
know that Hp^ i : LI' H is one-to-one and onto, thus is an isomorphism. It then 
follows that there exists a linear map F : H ^ Di such that for all z E H' , there 
exists unique x £ H such that 

(4.26) z = x + F{x) 

Claim 4.12. G{H') n H'^ = G{H') n 
Proof of Claim 4.12: For aU y € G{H') n H'^, z e iJ', we have y G G{H') C P™,;, 
thus 

(4.27) iy,z) = {y,x + Fix)) = {y,x)=0, 
which shows that y G LL^, so we obtain 

(4.28) G{H') n H'^ c G{H') n 
Now, for aU y' e G{H') n H^, x e H,we have 

(4.29) (y', x) = (y', x + Fix)) = (y' , z) = 
Thus 

(4.30) G{H')nH^ cG{H')r\H'^ 
Combining (|4.28p and (|4.30l) . we obtain 



(4.31) 



G{H') n iJ' - G(ff') n ff^ 



Q.E.D. 



steady state of chemical reaction networks 



17 



Now utilizing Claim 4.12 we obtain 

(4.32) G{H') n H'^ = G{H') n = G o {I + F){H) DH^^id 

where the last statement follows from Lemma 4.5. 
□ 

Lemma 4.13. If the CR has exactly I linkage classes, then for any subspace H of 
i?™, we have that 

(4.33) G{H) r\H^ ^% 

Proof. If the CR has I linkage classes, by relabeling Fi, I2, • • ■ , Ym we may assume 
that R has the same diagonal block form as that of (|2.2p . 

Now let H C i?" be fixed. By linear algebra, we know that H ^ H'-^^ © H'^^\ 
where H^^'> = H f] Di C Di,H^'^'> Ci Di = {0}. Thus we have 

\ =G(ffW®i7(2))ni/(2)^ 
(4.34) 

since combining H^^^ C Di = P,^ , and G(i?™) C Pm,i we have G{H^^'> ® H^^)) c 
Pjnd C H'-^^'^ . Finally noting that 

G(i7(2))ni/(2)^ ^0 

due to Lemma 4.9, we have 

G(i7(i) © ^ jj(2)^ 3 G(i/(2)) n ^ 0. 

□ 

5. Proof of Main Result and the index formula. Now we only need to 
observe that the existence problem for each compatibility class can be reduced to one 
of intersection problem type as in the previous section. 

Proof of Theorem 1.1: First a few words about the notation and definitions. 

Let G, R, i^T, be defined as in Section 2. Suppose that the CR has I linkage 
classes, by relabeling Yi, I2, ■ • • , Ym we may assume that the relation matrix R has 
the same diagonal block form as that of (|2.2p . Let S be the stoichiometric subspace 
in i?", and define W{xo) = ln((a;o + S) n i?!^) -C* (ZK, where xq € R^. 

As in the proof of Lemma 4.13, we have the decomposition of K as K = Ki © 
where Ki ^ K r\ Di a DuK'^C] Di ^ {0}. Now let K2 = Hp^ , {K'2). Then similar to 
the proof of Lemma 4.11 we have 



(5.1) 



G{W{xo)) r\K^ ^ G{W{xo)) n K'^^ ^ G{W{xo)) n K2 
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Next we will need the following lemma, for the proof, the reader is referred to 
Proposition B.l in (0]). 

Lemma 5.1. Given xq G i?", we have : (xq + S) i?" — >■ K2 given by 
^(a;) = 11/^2 (In a; • C*) is one-to-one and onto, actually, a dijjeomorphisni. 

Due to Lemma 5.1 and the fact that W{xo) C K, Hp^^^ = K2, there exists a 
continuous ( probably nonlinear) map F : K2 — > Di such that for all x € W{xo), 
there exists unique z G K2 such that x = z-\-F{z). Defining a function G : K2 — > Pm,i 
as 



(5.2) G{z) = G{z + F{z)) 

for aU z e K2. Note that G{K2) = G{W{xo)). 

Due to Remark 4.9 (or Corollary 4.10) of Section 4 we know that 



(5.3) G(i^2) n if2^ = G o (/ + F){K2) n K2^ ^ 0, 

i.e., we have G(W{xo)) n i^s^ 7^ 0- 

Combining ((5?T|) with ((Ot we have that G(W(xo)) n A'^ 7^ 0. Now we only need 
to observe that the map F we defined above is an analytic function, thus by Remark 
4.8 of Section 4, we must have car{G{W{xo)) n K-^) is finite, where car{Q) means 
the cardinality of the set Q. The proof of Theorem 1.1 is complete. 

Corollary 5.2. Let s be the dimension of the stiochiometric subspace. If we 
restrict the vector field f in U.l]) to some fixed positive stiochiometric compatibility 
class, denote the steady states in the stoichiometric compatibility class as xi,X2, ■■■ ,Xt 
(the number t of steady states depends on the stoichiometric compatibility class we 
choose), then we have 



J2ind{x,) ^ {-ly 



i=l 



where ind{xi) is the index of vector field f at steady state Xi,i = 1,2,. 
Proof. Let G, R, K, K2, S be defined as above. 
First, we give the following 



Definition 5.3. For the configuration matrix G — 
reduced configuration matrix G as 



, we define the 



G = G- 



G 
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where Di = ■ 



1 1 
1 1 



1\ 



, with rui given at the beginning of section 



V 1 1 ... 1 / 
4 for i = 1,2,...,Z. 

It is easy to check that C has the foUowing property 



(5.4) 



R C = R C 

nK,(x-c*) = nK,(x-c*) 



Now fix xq G i?" , from the proof of Theorem 1 . 1 above we know the vector field 
/ on xo + 5 can be written as 



(5.5) 



dx 



where z — "^{x) = H-k^ (In a; • C*) € K2. By Lemma 5.1 we know that ^ : + S* — > K2 
is a diffeomorphism. Thus considering the induced vector field ^I'* o / on K2 we have 



(5.6) 



f =vE,,o/oxI/-i(z) = nK,(f .Dza5(i)-C*) 



nif,(e-+^(^) •R-C-i?iag(i)-C*) 



where Diag{^) = 
reduced configuration matrix C. 



X2 



\ 



Now notice that the matrix C • 



V 



and we have used property (|5.4p of the 



j_ 

2:2 



•C* is positive definite when 



the action is restricted to i^2, thus by a simple homotopy between C- 



^ ^1 



V 



X2 



(7* and the identity matrix we see that the vector field \I>* o / o ^f"^ is homotopic to 
nif2(e^^^^^-' • R), while the critical points and corresponding indices remain invariant 
under this homotopy. 

We observe that in the proof of Lemma 4.5 we construct a bounded, convex 
domain Q so that the vector field ll^^ (e^+^*^^' • R) is pointing inwards at each point 
of dV,, thus by a standard result of degree theory we have that 



t t t 
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where s is the dimension of K2, which by Lemma 5.1 is equal to the dimension of S, 

the stoichiometric subspace. 

□ 
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